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We study two simple mesoscopic models of interacting two electrons; first one consists of two 
quantum coherent parallel conductors with long-range Coulomb interaction in some localized region 
and the other is of an interacting quantum dot (QD) side-coupled to a noninteracting quantum wire. 
We evaluate exact two-particle scattering matrix as well as two-particle current which are relevant 
in a two-particle scattering experiment in these models. Finally we show that the on-site repulsive 
interaction in the QD filters out the spin-singlet two-electron state from the mixed two-electron 
input states in the side-coupled QD model. 
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| /' | | The Landauer-Biittiker (LB) scattering approach is the cornerstone in the study of quantum transport in nonin- 
teracting mesoscopic systems [H, [2]. One can make an one-to-one connection between the Lippman-Schwinger (LS) 
scattering theory and the LB approach [3|. There are also several theoretical approaches to incorporate Coulomb 
interaction between electrons to investigate the transport phenomena in interacting models 0, But most 

i-C ' of these techniques are either perturbative in the interaction/tunneling strength or valid only in the linear response 
regime. Thus a full-fledged quantum transport method to study the interplay between the strong interaction and the 
noncquilibrium behavior is on-demand. One way to tackle this problem is to employ the time-independent elastic 
LS scattering theory and find an exact many-body scattering eigenstate of the open interacting system. The basic 
assumption here again as the original LB scattering approach is that all the dissipation is considered to occur only 
q-i in the reservoirs connected to the mesoscopic sample. Recently there have been some studies [H, [13, El, EH, 0, EH 
along this direction. A model of two quantum coherent conductors interacting weakly via a long range Coulomb force 
locally in some region has been studied in Both the LS and the LB approaches have been employed and the 

two-particle scattering matrix is expressed in terms of the scattering matrices of the noninteracting conductors. The 
results in [1 01 ] are perturbative in the interaction strength. Here we study a similar model and show that it is possible 
to find exactly the two-particle scattering matrix as well as the two-particle current change due to the interaction 
in this model. Later we investigate another interacting open quantum impurity model; an interacting quantum dot 
(QD) is side-coupled to a noninteracting quantum wire. We show that the on-site interaction in the QD filters out 
If} . the two-particle spin-singlet state from the mixed two-particle input states. We here apply the technique developed 
in Ref.fll], [II] based on the LS scattering theory. It has been shown in [ll], [T§] that one can find an exact two- 
0^ . particle scattering state for certain open quantum impurity models. A many-particle scattering state has been found 
in [ill . [HJ within a two-particle scattering approximation. Physically in a real two-particle scattering experiment [TBj 
, one considers two wavepackets representing two electrons. 

II. SCATTERING OF ELECTRONS BETWEEN TWO INTERACTING CONDUCTORS: 

We consider two quantum dots capacitively coupled via interaction V. Both the dots are connected to two nonin- 
teracting leads modeled by one-dimensional tight- binding Hamiltonian, 71%. Electron moves from one lead to other 
through the dot and interacts with electron of the other dot only at the dot sites. But there is no exchange of electrons 
between the dots. This is a lattice version of the model studied in [io| . We can better think of the model as of two 
separate parallel mesoscopic conductors (labelled by / and //) in proximity of each other and single electron in each 
conductor (see FigfT]). Electrons in the conductors interact only in some localized region. For simplicity we consider 
here spinless electrons. The Hamilton 

H = U 1 + H 11 + V, where U a = Hi + H a D + V a with (2.1) 
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FIG. 1: (color online) A schematic description of the two coherent parallel conductors in the presence of finite on-site energy 
at the dot sites. 



Above 2~D' implies omission of I = —1,0 from the summation and a = I, II. Here c a> i (ci ; ) is the electron annihilation 
(creation) operator in the ath conductor. e Q is the on-site energy on the ath QD and A is the strength of electrostatic 
Coulomb interaction between electrons in the two QDs. Also we set the hopping amplitude between the sites on both 
the conductors to unity. The double QDs with a purely capacitive interdot interaction can be labeled by a pseudospin 
index for the two dots and thus can be considered as a realization of the Anderson impurity model. One expects 
that electron transport through the QDs, that are weakly tunnel coupled to their leads, is dominated by the interdot 
interaction at low temperatures and this leads to Coulomb blockade jfj|. Recently Hiibel et al. have shown that the 
interdot Coulomb blockade can be overcome by correlated tunneling when tunnel coupling to the leads is increased 
We have here complete freedom within our approach to tune the tunnel junctions as well as the values of the 
dot energies and interdot interaction. Therefore we are able to study all regimes of the parameter space in our work. 



A. Scattering states 



We find exactly all the two-electron energy eigenstates for this model. First we need to calculate two-particle 
eigenstates of the noninteracting Hamiltonian Ha = Ti 1 + TL 11 for A = 0. One-electron eigenstates 4>%(l) (= {l\<f>%)) 
in a single conductor of Hamiltonian TL a (with a — I, II) can be found by solving the single electron Schrodinger 
equation. For an electron incoming from the left < k < 7r, the complete wave function is given by 

= e lkl + r%e~ m for I < -1, 
= (l + r£)/ 7 Q forZ = 0, 

= t%e iM for I > 1, (2.2) 
with the following transmission (reflection) amplitudes t% (r^), 

, a 2i7 Q2 sinfc 

2^ 6 — £q — 2 cos k fc fc 

Similarly one can find the single particle scattering state for a particle incident from the right. 

We form a two-particle incoming state 0k(l) of two electrons with one in each conductor. 0k(l) = (h)^ 1 (h), 
with k = (feijfe) and 1 = (h,h)- As electrons in the two conductors are distinguishable, we need not to anti- 
symmetrize the two-electron wave function. The energy of this state is = + E^ 2 . A scattering eigenstate |^k) 
of H with energy E^ is related to a state |</>k) (= |$L)|0fc J )) of "Ho through the Lippman-Schwinger equation Q 

l^> = I'M + C? ± (i? k )V|^}, (2.3) 
where G^(E^) = 



Eu — Hoi: it 

As usual, the indices (±) indicate outgoing wave (+) or incoming wave (— ) boundary conditions. Now in the two- 
electron sector, with the position basis |1) and an incident state < l\<f>) = </>k(l), Eq. (|2.3|) gives 



V+(l) = k (l) + XK Ek (l) ^+(0), (2.4) 
+ < 



where K Ek (l) = < l|G+(£ k )|fj> 
and = (0, 0). We can determine ip£ (0) using Eq. (|2.4|) . 



(2.5) 
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The two-electron scattering eigenstate is completely given by Eqs. (| 2 . 4ti 2 - 5 1) . As it has been shown in [TT|, [l2j that 
after scattering from the interaction the total momentum of the scatterred particles is not conserved though the total 
energy remains same. The momenta (k[,k' 2 ) of the scatterred particles are related to the incident momenta (ki,k 2 ) 
through, cosfci + cos k 2 = cosfc^ + cosk' 2 . The matrix elements Ke^(X) are known explicitly and are given by 

KEjl) = £ £ W ^ - k + t f 0^g W(°)C(0)- (2-6) 



B. Exact scattering matrix 

We define the two-particle scattering matrix following [l(|; here we suppress the lead index. The two particle- 
scattering matrix is given by 

S(E ll E 2 ,E 3 ,E 4 )S(E 1 +E 2 -E 3 -E 4 ) = (^eJ^eJ- ( 2 -7) 

The energy or momentum indices in the two-particle outgoing state indicate the energy or momentum of the two 
incident electrons. After some rearrangement using the Lippman-Schwinger equation, we find 

WeuEiWe^Ei) = (<l>EuE2 \<I>E 3 ,E 4 ) + (4>E U E 3 \ TgT + E^- H Q + it) V ^ 3 ^ 

+ (E 1 +E 2 -l 3 -E 4 + l e) { ^- m ^ ) 

= 8{E 1 - E 3 )S(E 2 - Ei) - 27Ti(^ 1 , E2 |V|^ 3 , E4 ). ( 2 - 8 ) 
Finally we calculate the change of the two-particle scattering matrix due to the interaction, 

5S(E 1 ,E 2 ,E 3 ,E 4 ) = -27Ti(^K^ 1 |V|^ 3 , E4 ) 

= -2^(^0)^10) (0,0^+^) 

1 - AA(£ 3i£ ; 4 )(U) 

where we have used Eq. (|2.5|) in the last line. In the weak coupling limit, i.e., A — > 0, one gets back the corresponding 
expression of the two-particle scattering matrix of [10]. Eq. (|2~9p is one main result of this paper. We emphasize that 
due to the interaction the two particles can exchange energy after scattering. 



C. Two-particle current 

The current density in the conductor I is given by the expectation value of the operator, jf = — i{c\ci + i — h.< 
the two-electron scattering state |^>k) 



1 0k) + 1 5k) (from Eq. (|2.3p ). The current in the incident state is given by 

(Mli\M = 2Afx\t I kl \ 2 sm(k 1 ), (2.10) 

where Ax (a normaliasation factor) is the total number of sites in the conductor I. Similarly, one can find the current 
in the conductor II in the case A = 0, (0k|ij ! 10k) = 2A/xx|^ 2 sin(fc2), where Afxx is the total number of sites in the 
conductor //. The change in the current in the conductor / due to scattering, 5j I (ki,k 2 ) = (^kb'/IV'k) — (0k|j/|0k), 
gets contributions from two parts, namely, jg = (S^jl \Su) and jjj = (S'k|j; / |0k) + (0k|i/|<Sk)- 



j| = 2 Im[<£ k |cJcj+i|Sk)] - 21m 



A 2 |^+(0)| 5 



2tt 



dqi I (qi)I*(qi) 



with I s ( qi ) = 1- f dq ^ (0) . 4£(l + 8) with s = 0,1 
2?r £ k - E qqi ~ie q 



(2.11) 
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FIG. 2: (color online) Plot of the interaction induced two-particle current change in conductor / with time (t) for different 
tunneling strength. Here fci = fe = 1.35, <j> = 0.35, A = 0.8 and id = 0.01. 



If we switch off the on-site energy in the two interacting dot sites and take the hopping energy indentical to that of 
the leads, i.e., ej = ejj = 0, "f 1 — J 11 = 1, then we are able to integrate the Eqs. (|2. 1 112. 12)) analytically and the total 
two-particle current change is given by 

Sj T (h,k 2 ) = ! jy^ftL [sgn(fei) +sgn(fc 2 )] (2.13) 

Thus far we could not calculate Eqs. (|2. 1 112. 12jl analytically for arbitrary values of ej, eij,^ 1 , j 11 , instead we evaluate 
them numerically [ll[ . We find that the two-particle current change due to the interaction is smaller by a factor Mi 
than the incident current; this signifies that the probability of two-electron collision in conductor / is order of 1/A/j. 



D. Periodically varying on-site energy 

In a recent experiment [TtJ the noise cross-correlation of two capacitively coupled QDs in the Coulomb blockade 
regime has been measured and the sign of this correlation has been found to change sign with tuning the on-site energy 
of the dot site by the gate voltage. Inspired by this experiment we now evaluate the two-particle current change for a 
periodically varying on-site energy of the two dots. This is like a two-electron quantum pump and we wish to study 
exactly how the strength of the interaction or the tunneling affect the current change averaged over a full cycle. We 
use, 

ei(t) = cos(wt) and, £u{t) — cos(uit + <fi) (2-14) 

In Fig[21 we plot the interaction induced two-particle current change 5j(ki,k 2 ) with time for different values of the 
tunneling between the dot and the leads. One can understand qualitatively that for which values of the on-site 
energy the sign of Sj (ki,k 2 ) is changing. In the absence of the interaction (A = 0), for a fixed incident energy 
Ek x (= — 2cosfci), a single particle resonance occurs at the on-site energy of e/, given by 

ei = E kl - y 2 (E kl T 05 - 4) (2.15) 

In the weak coupling limit J 1 — > we expect £/ = . Indeed we find that the change in sign of Sj (ki, k 2 ) (like in 
FigJ5]for the weak coupling case) occurs for the on-site energies correspond to the single particle resonance. 

Also we see in Fig[3]that the average two-particle current ({6j I (ki,k2))) depends on the interaction strength and 
can be positive or negative depending on the interaction. (Sj 1 (ki, k 2 )) crosses over from a positive to a negative value 
as the interaction strength increases. We define 



(5j I {k l ,k 2 )) = ^ j dtfik^m 



(2.16) 
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FIG. 3: (color online) Plot of the interaction induced two-particle current change in conductor / with time (t) for different 
interaction. Here ki = fe = 1.35, cf> = 0.35, ■y 1 = r y 11 = 0.6 and w = 0.01. 



Finally we write down formally the expressions (24-26) of Ref.[lfJ for arbitrary interaction. 

(ftjj) = (V'fc 1 ,fc 2 |4 c ilV'fc 1 ,fc 2 ) = ^2 \ip ku k 2 (l,m)\ 2 , 1 > 0, ki > 

m 

(^r^r) = ( ,( /'fc 1 ,fc 2 |c| c i rf L^m|V'fc 1 ,fc 2 ) = |^fc 1 ,fc 2 (Z,m)| 2 , l,m > 0, ki,k 2 > 

= " (2-17) 

As the two-particle scattering state, tpk 1 ,k 2 (^ m ) * s known explicitly from Eqs. (|2.4fl2.5|) . one can calculate these corre- 
lations for any A. 



III. A SIDE-COUPLED INTERACTING QUANTUM DOT ACTING AS TWO-ELECTRON SPIN 

FILTER 



In the parallel conductors model one expects that the two electrons in the different conductors get entangled 
(orbital or pseudospin entanglement) due to the interaction. Here we study another mesoscopic system with the 
localized interaction which acts as a two-electron spin filter, i.e., the side-coupled interacting quantum dot filters out 
a two-electron spin-singlet state in the output lead from two-electron mixed input states in the input lead of the 
noninteracting quantum wire. Recently there is one similar study with the Anderson impurity model for the linear 
energy-momentum dispersion of the leads [l8l |. 

We consider an interacting QD side-coupled to a perfect quantum wire (QW) [see FigU] modeled by a single 
electron tight-binding Hamiltonian. The dot consists of a single, spin-degenerate energy level with an on-site Coulomb 
interaction between electrons. The main idea of our scheme is to prevent single-electron tunneling as well as current in 
the spin-triplet channel in the output lead. In that sense our program here matches with that of Oliver et al. fTi) . But 
we achieve these criteria through different mechanism and here, we don't need a three-port quantum dot geometry as 
well as leads acting as "energy filters". We avert single-electron tunneling in the output lead by tuning the voltage 
gate attached to the QD. It occurs due to destructive interference of the single electron wave when the on-site energy 
of the dot-site is same as the energy of the incident electrons. It turns out that the above condition is also sufficient 
for complete destructive interference in the spin-triplet state. On the other hand in the presence of a finite Coulomb 
interaction in the dot, we get a finite current solely comprises of the two-electron spin-singlet state. Thus exchange 
interaction and quantum interference mediate to filter out the spin-singlet state of a two-electron mixed-state input 
of opposite spins. The full Hamiltonian TL of the system consists of a noninteracting part TLq and an on-site Coulomb 
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FIG. 4: (color online) A schematic description of a quantum wire with a side-coupled quantum dot (QD) modeled as an 
Anderson impurity. 

interaction part V. 

H = Ho + V 

OO 

(= — 00,CT=t,i f=Tll 
<T=T,J. 

V = Un d] n di . (3.1) 

where n^.o- = c<j )CT is the number operator in the dot for spin a and e d is the on-site dot energy. U is the strength 
of the on-site Coulomb energy in the dot site and 7 represents the tunneling strength between the quantum wire and 
the quantum dot. Again we set the hopping amplitude between sites on the quantum wire to unity. 



A. Scattering states 

For an electron coming from the left, the eigenstates \4>1) of Hq in the position basis \l) are given by 

4,1(1) = e lkl + r k e~ lkl for / < 0, 

= ^-^ on the dot site, 

ed + 2 cos k 

= t k e iM for I > 1, (3.2) 

where <pZ(l) = !</>£) an( i < k < n (with a =|, j). The transmission, reflection amplitudes t k , r k are determined by 
solving the single electron Schrodinger equation; they are, 



2i(2 cos k + ed) sin k 
2i(2 cos k + ed) sin k + j 2 



The incident energy of a single electron is E k — — 2cosfc. From the Eqs. (|3.3p , we see that for a finite 7, the 
transmission amplitude vanishes at a finite incident energy, E k = td, i-e., when the on-site dot energy is same as 
the energy of the incident electron. We can achieve this criterion by controlling the plunger gate acting on the 
QD. As before we calculate the single electron tunneling current by taking expectation value of the current operator 
ji = — i{c, f _c;_|_x,<r — h.c.) in the single electron scattering state \4>D- Then, the single electron tunneling current in 
the output lead is given by 2A/"|££| 2 sinfc, which vanishes at E k = £d- Here TV is a normalisation factor depicting the 
total number of sites in the entire system. Now we consider that a spin-up and a spin-down electron incident in the 
input lead of the QW. The two-electron input state (with total S z — 0) is a mixture of spin-singlet and spin-triplet 
states whose spatial wave-functions are respectively symmetric and anti-symmetric. So the on-site Coulomb repulsion 
can cause scattering between two electrons in the spin-singlet channel but not in the spin-triplet channel. We find the 
current contribution from the spin-triplet channel in the output lead by taking expectation of ~^2 a 3° in the spin-triplet 
scattering state. If momentum of the two incident electrons are (^1,^2), the current in the spin-triplet channel is 
given by, 2N(\t\ i | 2 sin(fci) + \t k2 1 2 sin(fc 2 ))- Then for vanishing current in the triplet channel we need to satisfy, 
E kl = Ek 2 = f-d- The last criterion also eliminates the possibility of both up or down spins (with total S z — ±1) in 
the spin-triplet part of the input channel. 

Now we calculate the effect of Coulomb interaction on the scattering of electrons in the spin-singlet channel. As 
there is no spin-flip interaction in the Ti., we need to consider only the spatial part of the spin-singlet wave function. 



7 



^ 1 ■so 


/ 

/ 

/ 

/ 




C 

1) 
p 

u ioo 


/ 

/ 
/ 


-- y=0.3 

— y — u.j 

— 7=0.8 








"m 
_« 


l 




S 50 


1 




& 

GO 


1 X 
( 





Interaction (U) 



FIG. 5: (color online) Plot of the spin-singlet current arising from Js = ^ CT j% (feo, ko)/2 vs interaction ({/) for = —0.2 and 
different values of the coupling strength 7. Inset shows a plot of j° (ko, fco) with U for three different values of td and a fixed 
7 = 0.3. In the inset, the upper curve is for td = —0.2 and the lower two (almost overlapped) curves are for ed = —0.6, 0.6 
respectively. 



The scattering of two electrons in the spin-singlet channel can be studied using the LS formalism of Ref . [llj . We 
consider here coherent electron transport at zero temperature. Incoming state of two electrons in the position basis is 
given by (fk(l) = <p' k {h)<j> k2 (h) with k = (fei, fe) and 1 = (ii, Za). Then using the LS equation we find the two-electron 
scattering eigenstate \ip) of H as 

tfk(l) = ^ k (l) + UK Ek (\) Vk(d), (3.4) 
where K Ek (l) = (l|G+(£ k )|d> and ^(d) = • 

with d = (d, d) . Again the two-electron scattering eigenstate of Ti. is completely given by Eq. (|3.4p with 

^(i) = f r p — L— -^ q (i)<(d). (3.5) 



B. Spin-singlet current 



When the energy of the incident electrons is same as the on-site energy of the QD, there is no single electron tunneling 
as well as a zero current in the spin-triplet channel. So the current in the output lead is solely determined by the 
contribution from the spin-singlet channel. Now we determine the two-electron current density by the expectation 
value of the operator 53 CT if in the scattering state ?/>k) = ¥>k) + 5k) (from Eq. (|3.4p ). We calculate different parts of 
the current in the spin-singlet channel separately. The current in the incident state (with incident wave- vector ki,k 2 ) 
is given by 



(¥>k|$>>k) = 2A/-(|4j 2 sin(fc 1 ) + |4j 2 sin(fc 2 )), 



(3.6) 



which vanishes for E kl = E k2 = td identically. The change in the current due to the interaction, dj a (ki, k 2 ) = 



(V'kbTlV'k) - <<Pk|jf>k) = j° + f c , with f s = (S k \j[\S k ) and ft = {S k \tf\y k ) 
the current change solely from the scattered wave-function |Sk) is given by 



{ipk\ii\Sk) ■ For electron with "f spin, 



(SkljhSy) = 2Im[[/> k (d 



dq 



dg 2 <(Z + l)^*(ri ) fdq 3 4>ll{l)4>lM 
2ir E^-E qiq2 +i 



2tt E u 



E qiq3 it 



(3.7) 



This expression is similar in form of Eq. (|2.1ip . Now ^k(d) is nonzero even for fci = k 2 = cos x (— e^/2) = kg and the 
integrand of Eq. (|3.7p can not be said to be zero a priori. So we expect to have a finite contribution in 5j a {ko, fco) from 
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this part. The other term in the current change comes from the overlap between \Sk) and |vk), which is given by 

The factors <Pk*Q) an d Vk*^ + -0 m Eq. (|3.8[) vanish for fei = fco hi the output lead, i.e., / > 0. So there is no 
contribution in Sj a (kg, ko) from the term in Eq. (|3.8p if we evaluate current change in the output lead for electrons 
being incident from I < 0. Ultimately we need to evaluate the integral in Eq. (|3.7p to quantify the amount of spin- 
singlet pair generated in the output lead. As the parallel conductors model we determine it numerically for different 
values of the coupling strength 7 and the on-site dot energy e^. We plot Js(fco, fco) = Y]^ j% (fco, fco)/2 with interaction 
U for different 7 in Fig[5l which shows that the spin-singlet current increases with weaker coupling of the QD with 
the transport channel. This can be understood from the Eq. (|3.7[) . jg(k ,k ) depends on | -0k o | 2 which is inversely 
proportinal to j 4 . Occupation probability of the singlet pair at the QD is higher for smaller 7; so the electrons 
scatter strongly with smaller 7 and (fco, fco) increases. FigJSjalso shows that Js(fco, fco) saturates after some critical 
strength of interaction, U c which becomes smaller with decreasing 7. Here we should also clarify that one needs a 
finite coupling of the QD with the quantum wire to get a antiresonance in the single electron tunneling. We plot 
Js(fco, fco) for three different values of ed in the inset of FigO we find that the magnitude of Js(fco, fco) is same for a 
dot and an antidot on-site energy and the current increases with smaller dot energy. To check that the total current 
is same after scattering in both the input and the output leads, we evaluate <5j CT (fco,fco) hi the input lead also, i.e., 
I < 0, in which case we need to evaluate both Eqs. (|3.7l3.8p . We find total current is same in the input and the output 
leads within small numerical error. 

IV. CONCLUSION 

We have calculated exactly the two-particle scattering state as well as the corresponding current in two interact- 
ing mesoscopic lattice models. In principle one needs to find a many-particle scattering state to study the out of 
equilibrium phenomena in these impurity models. But recently it has been shown in [Til , [lil ll8T ] that many of the 
nonequilibrium quantities like the current-voltage characteristics have significant features in the two-particle current 
for weak interaction or low density of electrons. Though the many body effect drastically changes for strong interaction 
or higher density, for example, one expects to find an anti-Kondo resonance in the conductance of the side-coupled 
dot model in the presence of many electrons in the quantum wire. 

We thank D. Sen and M. Buttiker for many fruitful discussions and useful suggestions on the draft. The hospitality 
of Dept. of Theoretical Physics, Unversity of Geneva is gratefully acknowledged. 
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